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Abstract 

This is an introductory document surveying several results in polynomial approximation, 
known as the Bramble-Hilbert lemma. 

1 Introduction 

In numerical analysis, the Bramble-Hilbert lemma bounds the error of an approximation of a 
function u by a polynomial of order at most m — 1 in terms of derivatives of u of order m. Both the 
error of the approximation and the derivatives of u are measured by LP norms on a bounded domain 
in W 1 . This is similar to classical numerical analysis, where, for example, the error of interpolation 
u on an interval by a linear function (that is, approximation by a polynomial of order one) can 
be bounded using the second derivative of u. The difference is that the Bramble-Hilbert lemma 
applies in any number of dimensions, not just one dimension, and the approximation error and the 
derivatives of u are measured by more general norms involving averages, not just the maximum 
norm. 

Additional assumptions on the domain are needed for the Bramble-Hilbert lemma to hold. 
Essentially, the boundary of the domain must be "reasonable". For example, domains that have 
a spike or a slit with zero angle at the tip are excluded. Domains that are reasonable enough 
include all convex domains and Lipschitz domains, which includes all domains with a continuously 
differentiable boundary. 

The main use of the Bramble-Hilbert lemma is to prove bounds on the error of interpolation of 
function u by an operator that preserves polynomials of order up to m— 1, in terms of the derivatives 
of u of order m. This is an essential step in error estimates for the finite element method. The 
Bramble-Hilbert lemma is applied there on the domain consisting of one element. 



2 The one dimensional case 



Before stating the lemma in full generality, it is useful to look at some simple special cases. In one 
dimension and for a function u that has m derivatives on interval (a, b), the lemma reduces to 

inf \\ u ( k ) - v ( k )\\ < C{m){b-a) m - k \\u^\\ JP( M , 

w6 p m _ l ll \\LP(a,b) - ^ ' ^ I N ULP(a,b)> 

where P m -i is the space of all polynomials of order at most m — 1. 
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In the case when p = oo, m = 2, k = 1, and u is twice differentiable, this means that there 
exists a polynomial v of degree one such that for all x € (a,b), 

\u (x) — v (x)\ <C{b — a) 2 sup \ u"\ , for all x S (a, 6) . 

(a,6) 

This inequality follows from the well-known error estimate for linear interpolation by choosing v as 
the linear interpolant of u. 



3 Statement of the lemma 

Suppose O is a bounded domain in R n , n > 1, with boundary and diameter d. Wp(Cl) is the 
Sobolev space of all function u on Q with weak derivatives D a u of order |a| up to k in L p (f2). 

Here, a = (ati, ai, ■ ■ ■ , a n ) is a multiindex, \at\ = a\ + 02 ^ 1- a ra , and D Q denotes the derivative 

a\ times with respect to x%, «2 times with respect to «2> and so on. The Sobolev seminorm on 
W™(0,) consists of the L p norms of the highest order derivatives, 

i/p 

Ww^in) = I Yl \\ Dau Wlp(U) I if 1 < P < oo 

V |a|=m 



and 



HQ) — , max P a?i llL°°(Ji) 
jy ' \a\=m v ' 



Pk is the space of all polynomials of order up to k on W 1 . Note that D a v = for all v G P m -\- 
and |a| = m, so |n + v| Wm / n j has the same value for any v € Pk-i- 

Lemma (Bramble and Hilbert) Under additional assumptions on the domain 0, specified 
below, there exists a constant C = C (m, Q) independent of p and u such that for any u £ Wp(Q) 
there exists a polynomial v € P m -i such that for all A; = 0, . . . , m, 

\ u ~ v \\Vp(Q) - Cd m k \u\wp(n) ■ 



4 The original result 

The lemma was proved by Bramble and Hilbert [T] under the assumption that f2 satisfies the strong 
cone property; that is, there exists a finite open covering {Oi} of d£l and corresponding cones {Cj} 
with vertices at the origin such that x + C% is contained in f2 for any x G D 0{. 

The statement of the lemma here is a simple rewriting of the right-hand inequality stated in 
Theorem 1 in \X\. The actual statement in pQ is that the norm of the factorspace W™(Q)/P m -i 
is equivalent to the W™{Vt) seminorm. The W™{£1) is not the usual one but the terms are scaled 
with d so that the right-hand inequality in the equivalence of the seminorms comes out exactly as 
in the statement here. 

In the original result, the choice of the polynomial is not specified, and the dependence of the 
constant on the domain Q is not given either. 



5 A constructive form 

An alternative result was given by Dupont and Scott [1] under the assumption that the domain 
SI is star-shaped; that is, there exists a ball B such that for any x G fi, the closed convex hull of 
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{x} U B is a subset of fi. Suppose that p max is the supremum of the diameters of such balls. The 
ratio 7 = d/p max is called the chunkiness of fi. 

Given a fixed ball B as above, and a function u, the averaged Taylor polynomial Q' m u is defined 

as 

Q m u = J T™u(x)iP(y)dx, 

B 

where 

m— 1 1 

T-u (*) = ^ £ —D ol u (y) (x - y) a 

k=0 \a\=k 

is the Taylor polynomial of degree at most m — 1 of u centered at y evaluated at x, and ip > is a 
function that has derivatives of all orders, equals to zero outside of B, and such that 

ipdx = 1. 



Such function ip always exists. 

Then the lemma holds with the constant C = C (m, n, 7), that is, the constant depends on the 
domain f2 only through its chunkiness 7 and the dimension of the space n. For more details and a 
tutorial treatment, see the monograph by Brenner and Scott [2]. The result can be extended to the 
case when the domain £1 is the union of a finite number of star-shaped domains and more general 
polynomial spaces than the space of all polynomials up to a given degree [3]. 



6 Bound on linear functionals 

This result follows immediately from the above lemma, and it is also called sometimes the Bramble- 
Hilbert lemma, for example by Ciarlet [3]- It is essentially Theorem 2 from pp. 

Lemma Suppose that £ is a continuous linear functional on W™(0,) and ||^||^ym(Q)' its dual 
norm. Suppose that £ (v) = for all v E P m -i- Then there exists a constant C = C {Q) such that 

\£(u)\ < C \\£\\ W m(p,y \ u \w™{n) ■ 
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